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Effective Hamiltonians for LaFeAsO and LaFePO are derived from the downfolding scheme 
based on first-principles calculations and provide insights for newly discovered superconduc- 
tors in the family of LnFeAsOi-zFz, Ln = La, Ce, Pr, Nd, Sm, and Gd. Extended Hubbard 
Hamiltonians for five maximally localized Wannier orbitals per Fe are constructed dominantly 
from five-fold degenerate iron-3d bands. They contain parameters for effective Coulomb and 
exchange interactions screened by the polarization of other electrons away from the Fermi level. 
The onsite Coulomb interaction estimated as 2.2-3.3 eV is compared with the transfer integrals 
between the nearest-neighbor Fe-3d Wannier orbitals, 0.2-0.3 eV, indicating moderately strong 
electron correlation. The Hund's rule coupling is found to be 0.3-0.6 eV. The derived model 
offers a firm basis for further studies on physics of this family of materials. The effective mod- 
els for As and P compounds turn out to have very similar screened interactions with slightly 
narrower bandwidth for the As compound. 

KEYWORDS: first-principles calculation, effective Hamiltonian, downfolding, constrained RPA method, 
LaFeAsO, LaFePO, Oxypnictide, high-temperature superconductivity 



Recent discovery of a new superconductor, 
LaFeAs(OjF), 1 ^ has triggered extensive studies on 
the family of layered iron arsenide compounds with 
ZrCuSiAs-type or ThCr2Si2-type structure, whose 
superconducting critical temperature T c is now raising 
up to ^55 K. 2 ) A mother compound LaFeAsO shows 
antiferromagnetic order 3 ) with bad metallic transport 
properties 1 -* supporting a significant role of electron 
correlations. The antiferromagnetic ordered moment 
^0.36 hb is unexpectedly small, implying large quantum 
fluctuations arising from electron correlations with com- 
peting ground states. Superconductivity appears when 
carriers are doped by the substitution of F for O, 1 -* or 
the introduction of O vacancies. 4 > Many of experimental 
data support a significant role of electron correlation 
in realizing the superconductivity. 5 ~ 8 ^ Although the 
As compounds have the maximum T c ~ 56 K, the P 
compounds show one order of magnitude lower T c . 9 ) 

Conventional density-functional calculations with the 
local density approximation (LDA) or the generalized 
gradient approximation (GGA) have clarified entangled 
ten-band structure near the Fermi level mainly originat- 
ing from five-fold degenerate iron-3c? orbitals contained 
in each of two iron atoms in the unit cell. 10-14 ' For the 
mother material, the initial LDA calculation 10 ) predicted 
the nonmagnetic ground state in close proximity to a fer- 
romagnetic metal while the recent results show antifer- 
romagnetic order. n ~ 13 ) In particular, the stripe-type an- 
tiferromagnetic order is correctly reproduced. 12 ' 13 ) The 
calculated ordered moment obtained so far ranges be- 
tween 1.2 and 2.6 ^b, 1113, 15 ' in contrast to the tiny 
ordered moment discussed above. Broad peak structures 
of magnetic Lindhard function calculated by using the 
LDA/GGA Fermi surface suggest competitions of sev- 



eral different ordering tendencies. 14, 16 ~ 18 ) 

Although overall experimental results suggest notice- 
able correlation effects, realistic roles of electron correla- 
tions on theoretical grounds are not clear. The relevance 
of the correlation effect is under active debates. 19 ~ 21 ) It is 
imperative to estimate at least the effective Coulomb in- 
teraction from first principles. Furthermore, since all the 
3d bands of Fe are as a first look wholly involved near 
the Fermi level, it is important to elucidate interplay of 
orbital degeneracy and electron correlation, which can be 
studied only by a model for degenerate bands. A reliable 
theoretical model derived for this family of compounds 
thus provides us with a firm basis for understanding the 
superconductivity and also a starting point for exploring 
further possibility of higher T c compounds. 

In this letter, we present ab initio low-energy effective 
Hamiltonians of LaFeAsO and LaFePO. Implications for 
the superconductivity in F-doped materials are also dis- 
cussed. A reliable downfolding scheme has recently been 
established which has enabled derivation of low-energy 
effective Hamiltonians from ab initio density-functional 
calculation of real materials. 22 ~ 24 ) The low-energy Hamil- 
tonian is derived after eliminating higher-energy degrees 
of freedom and estimating the renormalization effect 
from the high-energy part onto the bands near the Fermi 
level. The accuracy of the downfolding procedure has 
been established in various cases. 24 ~ 26 ) Here, we employ 
this downfolding scheme in deriving the extended Hub- 
bard models of LaFeAsO and LaFePO, consisting of band 
dispersion (kinetic energy) of electrons at maximally lo- 
calized Wannier orbitals (MLWOs), 27 ) as well as screened 
Coulomb and exchange interactions. Since the ten-fold 
Fe-3cZ bands are basically isolated from other bands, we 
derive effective models for these ten bands. 
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The extended ten-band Hubbard model reads 

U = ^ ^2 E *™R.nR.' a ™R a mR' 



a FLFL' nm 



+ 9 E E E] U mB.nB.'anR a mR' a mR' a nR 



<yp RR' nm 

+ JmRnR' (K^FpmR^nFpmR' + a n^R a nR a mR ia mR')^ > (!) 

where a° R (a^ R ) is a creation (annihilation) opera- 
tor of an electron with spin a in the nth MLWO cen- 
tered on Fe atoms in the unitcell at R. i m R„R' con- 
tains single-particle levels and transfer integrals, given 
by t mRnR , = {4>mR\7<o\4>nR<) with \</> nR ) = (J nR \0) and 
Ha being the one-body part of H. J7 m R„R' and J m R n R' 
are screened Coulomb and exchange interactions, respec- 
tively, expressed as 



t^mRnR' = {4>mR<t>mR | W\ 4> n R' </>„R' ) 



and 



J, 



mRnR' 



l R^nR'|W / |</>„R'0 m R) 



(2) 



(3) 



with W being a screened Coulomb interaction. There are 
already attempts to estimate i m R„R' using MLWO. 11, 14 ^ 
We here focus on an ab initio derivation of the many- 
body part of H; we estimate the interaction parameters 
in eqs. (2) and (3), based on a constrained random-phase 
approximation (cRPA). 22 ' 23 ) The cRPA has several ad- 
vantages over other methods such as constrained LDA. 28 ) 
We can precisely exclude screening processes among the 
Fe-3d MLWOs being the bases of the effective model. 
(This screening should be considered when we solve the 
effective model.) In addition, we can calculate matrix ele- 
ments of W as a function of R and R'; 29 ) i.e., we can ob- 
tain onsite and offsite interactions at one time. The cRPA 
becomes a good approximation, when the high-energy 
eliminated bands are well separated from the Fermi level 
and the screened Coulomb interactions between the high- 
energy electrons and the low-energy electrons are weak. 
As is noted above, quantitative accuracies of our down- 
folding including cRPA has already been confirmed in a 
number of examples. 22 ~ 25 > 29 ) In the present case of the 
iron compounds, the condition above is equally satisfied. 

U m RnR' of cq. (2) is practically calculated in the re- 
ciprocal space by using a Fourier transform of W, 

W{r,r')= J2 e j(q+G)r T^GG'(q)e^ (q+G ' )r '. (4) 

qGG' 

Here, G is a reciprocal lattice vector and q is a wave 
vector in the first Brillouin zone. We define W GG > (<z) as 



W GG ,(q) = 



4tt 



n \q + G\ 



GG 



q + G' 



where f2 is the crystal volume and e GG >{q) is the inverse 
dielectric matrix which is related to the irreducible po- 
larizability \ b y ecG'O?) = fee - v(q + G)xgg'(q), 
where v(q) = 47r/f2|q| 2 is the bare Coulomb interaction. 



The polarization matrix xgc (q) is calculated as 

X GG'(g)=EE'^«k +q | e i(q+G)r |# k ) 

k a/3 



-j(q+G)r 



l^ak+q} 



/ak+q ~ //3k 



where {ip a w} are the Bloch states and the prime attached 
to the band sum indicates that we exclude the 3d-3d band 
transitions in the calculation of \. By inserting eq. (4) 
into eq. (2), we obtain the form of 

A'jr 

U m0nR =- E e -^ R p m q(G)eaG'(g)p; q (G'),(5) 



qGG' 



where 



Pnq(G) 



N\q + G 



1 N 



; k+ q |e 



i(q+G)r 



IV'nk) 



with IV'nk) = X)r \4>nR)e lkR - It should be noted here 
that the quantities ( - 0mk+q|e i ^ q+G '' r |'i/'nk) can be easily 
evaluated with the fast Fourier transformation technique. 
Matrix elements of the bare (or unscreened) Coulomb 
interaction as t/^onR = {4>mO<i>mo\v\<l>nR(t>nB.) are cal- 
culated with replacing e GG ,(q) of eq. (5) by 6 GG >. The 
parallel treatment is applied to the derivation of screened 
exchange interactions in eq. (3). The result is 

JmOnR = -fr E P"^Rci(G)e GG ,(q)p^ nnRci (G') (6) 



qGG' 



with 



PmnRq(G) — 



N\q + G 



N 

E< 



ikR 



(^mk+q|e 



j(q+G)r 



|^nk>- 



We implemented this scheme in Tokyo Ab initio 
Program Package. 30 ^ With this program, electronic- 
structure calculations with the GGA exchange- 
correlation functional 31 ^ were performed using a 
plane-wave basis set and the Troullier-Martins norm- 
conserving pseudopotentials 32 ) in the Klcinman- 
Bylander representation. 33 ) Our iron pscudopotential 
was constructed under the reference configuration 
(3d) 60 (4s) 1 - 8 (4p)°° by employing the cutoff radii for 
the 3d, 4s, and 4p states at 2.1 Bohr. The energy cutoff 
was set to 64 Ry and a 9x9x5 fc-point sampling was 
employed. The experimental crystal-structure data were 
taken from ref. 3 for LaFeAsO and ref. 9 for LaFePO. 
The dielectric matrices were expanded in plane waves 
with an energy cutoff of 64 Ry and the total number of 
bands included in the sum in eq. (5) was set to 70. The 
sum of k in eq. (5) was evaluated by the tetrahedron 
method. The additional terms in the long-wavelength 
dielectric matrix due to nonlocal terms in the pseudopo- 
tentials was explicitly considered following ref. 34. The 
point q = G = in eqs. (5) and (6) requires special 
handling because of the singularity in the Coulomb 
interaction, which was treated in the manner described 
in ref. 34. 

The upper two panels of Fig. 1 shows ab initio band 
structures of LaFeAsO (left) and LaFePO (right). We see 
a good agreement between the original GGA band (red 
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line) and the Wannier- interpolated band (blue dots). The 
lower panels visualize our calculated MLWOs with the 
yz (left) and z 2 (right) symmetry. In principle, the band 
dispersion thus obtained should be renormalized by the 
interaction between electrons in the Fe-3(i and eliminated 
bands in the downfolding procedure. 22 ) In this letter, we 
assume that this self-energy effect is small and thus we 
employ the same dispersion for the downfolded Hamil- 
tonian. The overall band structure of LaFeAsO shows a 
feature similar to that of LaFePO with slightly (~ 20%) 
narrower bandwidth for the As compound. Around the 
r point, in addition to two hole cylinders with the d yz 
and d zx characters, a 2D-like d x 2_ y 2 band crosses near 
the Fermi level for the As compound, in contrast to a 
3D-like ribband of the P compound. 




Fig. 1. Upper panels: Ab initio band structures of LaFeAsO (left) 
and LaFePO (right). Red line and Blue dots are original-GGA 
and Wannier-interpolated bands, respectively. The zero of energy 
is the Fermi level. Lower panels: Isosurface contours of yz- (left) 
and z 2 - (right) MLWOs in LaFeAsO. The amplitudes of the con- 
tour surface are +1.5/-\A> (blue) and —1.5/^/v (red), where v is 
the volume of the primitive cell. Fe and As nuclei are illustrated 
by yellow and silver spheres, respectively. 



In Fig. 2, we plot diagonal elements of the macroscopic 
cRPA dielectric matrix, ej\i(q + G) = l/e GG (q), calcu- 
lated for LaFeAsO, as a function of \q + G\. We note 
that the resulting dielectric constant e° M at q + G — > 
exhibits a rather high value of 6.3, which is compared to 
that of transition metal oxides such as SrVC>3 (6.5). It is 
also interesting to note that €.m{q + G) does not severely 
depend on the direction of q + G, which means that the 
screening response is quite isotropic. 

Fig. 3 plots matrix elements of the screened Coulomb 
interaction, C/ m onR m eq. (5), denoted by green dots, for 
LaFeAsO as a function of the distance between the cen- 
ters of the MLWOs; r = \{<f>nR.\r\<f>nB.) ~ (4>mo\r\^>mo)\- 
In this plot, we set m to a d xy MLWO and display 
only the interactions between it and other MLWOs. 
For comparison, we also plot bare Coulomb interactions 




l 1.5 2 
\q+G\ (a.u.) 

Fig. 2. Macroscopic dielectric function of LaFeAsO as a function 
of \q + G\ obtained by the constrained RPA method. 



^mOnH as re d dots, which decay as 1/r (solid line) be- 
yond the nearest-neighbor Fe-Fe distance (>2.65 A). It is 
clear that the bare Coulomb interaction is significantly 
screened. In addition, as expected from Fig. 2, U m Q n ji 
decays as an isotropic function of l/(e° M r) (dotted line). 
Since the offsite interactions (E/mOnR^o) are more than 
five times smaller than the onsite U values, we may ne- 
glect them in the first step and just start with the onsite 
Hubbard model. The screened exchange interactions of 
JmOnTL in eq. (6) are found to decay very quickly; the 
magnitude is nearly zero, except for the onsite values. 




4 8 12 

r (Angstrom) 



Fig. 3. Calculated screened Coulomb interactions of LaFeAsO as 
a function of the distance between the centers of MLWOs. Only 
the interactions between a d xy MLWO at the homecell and other 
MLWOs are plotted. The red and green dots represent the bare 
and screened interactions, respectively. The solid and dotted 
curves denotes 1/r and l/(e9,r), respectively. 



We summarize in Table I the list of the important 
parameters in the low-energy effective Hamiltonian of 
eq. (1). The transfer integral between nearest-neighbor 
iron Wannier orbitals, t, is typically 0.2-0.3 eV, whereas 
the onsite screened Coulomb interaction, U, exhibits 2.2- 
3.3 eV. The strong orbital dependence comes from the 
fact that each orbital has different amount of leakage 
on the neighboring As atoms. The onsite exchange in- 
teraction (Hund's rule coupling) is 0.3-0.6 eV. We note 
that the resulting U and J values are smaller than the 
bare values by the factors of about 1/5 and 4/5, respec- 
tively. We note in passing that, for the two-degenerate 
dy Z and d zx orbitals, our computed J yz ,zx (= 0.45 eV) 



is close to the value estimated from (U y 



U y z,zx) /2 



(0.50 eV). The interaction parameters for the P com- 
pound are found to be very similar; for example, the 
onsite U values range from 1.9 eV for d x i_ y i to 3.3 eV 
for e? z 2 and d xy as in the case of the As compounds. 
The present value of U is substantially smaller than the 
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Table I. List of important parameters (in eV) in the present 
extended Hubbard Hamiltonian in cq. (1). From the top, the 
three 5x5 matrices represent transfer integrals between nearest- 
neighbor sites, onsite screened Coulomb interactions, and onsite 
exchange interactions of LaFeAsO. The bottom 5x5 matrix is 
transfer integrals of LaFcPO. The onsite energies for the five 
orbitals of LaFeAsO are (e X y, e yz , e z 2, e Z x, e x 2_ y 2) = (—0.12, 
+0.10, -0.19, +0.10, +0.16) cV. 



t (LaFeAsO) 


xy 


yz 


z 2 


zx 


x 2 -y 2 


xy 


-0.32 


-0.25 


-0.30 


-0.25 


0.00 


yz 


-0.25 


-0.21 


-0.08 


-0.13 


0.18 


z 2 


-0.30 


-0.08 


0.08 


-0.08 


0.00 


zx 


-0.25 


-0.13 


-0.08 


-0.21 


-0.18 


x 2 — y 2 


0.00 


0.18 


0.00 


-0.18 


-0.18 




[/(LaFeAsO) 


xy 


yz 


z 2 


zx 


x 2 — y 2 


xy 


3.31 


1.95 


1.89 


1.95 


2.09 


yz 


1.95 


2.77 


2.20 


1.78 


1.67 


z 2 


1.89 


2.20 


3.27 


2.20 


1.65 


zx 


1.95 


1.78 


2.20 


2.77 


1.67 


x 2 — y 2 


2.09 


1.67 


1.65 


1.67 


2.20 




J (LaFeAsO) 


xy 


yz 




zx 


x 2 -y 2 


xy 




0.54 


0.64 


0.54 


0.27 


yz 


0.54 




0.41 


0.45 


0.43 


z 2 


0.64 


0.41 




0.41 


0.50 


zx 


0.54 


0.45 


0.41 




0.43 


x 2 -y 2 


0.27 


0.43 


0.50 


0.43 






t (LaFcPO) 


xy 


yz 


z 2 


zx 


x 2 — y 2 


xy - 


-0.35 - 


-0.31 - 


-0.32 


-0.31 


0.00 


yz - 


-0.31 - 


-0.24 - 


-0.04 


-0.13 


0.18 


z 2 - 


-0.32 - 


-0.13 


0.13 


-0.04 


0.00 


zx - 


-0.31 - 


-0.13 - 


-0.04 


-0.24 


-0.18 


x 2 — y 2 


0.00 


0.18 


0.00 


-0.18 


-0.27 



values (~ 4 eV or larger) employed in model studies, 19 ' 20 ) 
which, as a first look, suggests that much smaller corre- 
lation effects should be expected in reality. However, it 
should be cautioned that the DMFT employed in the lit- 
erature 19 ' 20 ) ignoring spatial correlations may largely un- 
derestimate correlation effects. By considering the ratio 
U/t ~ 10 (the value itself is comparable to the cuprates) 
and the presence of the five-degenerate orbitals per Fe 
site, electron correlations in these compounds arc mod- 
erately strong. In addition, the present U values have 
orbital dependence with the order of 1 eV, which will 
generate nonnegligible orbital dependence of renormal- 
ization factor; it can be a critical parameter that affects 
the Fermi surface in quasipariticle band structure. Our 
downfoldcd effective Hamiltonian thus poses various con- 
straints on modeling of this family of materials and care- 
ful analyses by accurate low-energy solvers should be re- 
quired for reliable discussions of magnetic and supercon- 
ducting mechanisms. 

The effective models for LaFeAsO and LaFePO are 
basically similar in the band dispersion as well as in the 
screened Coulomb interaction. A main differences is the 
band width (see Fig.l) and this point can also be con- 
firmed from t values of LaFePO somewhat larger than 
those of LaFeAsO (see Table I). Another difference is a 
band close to the Fermi surface at T point; i.e., whether it 
has a two-dimensional dispersion with the d x i_ y i charac- 
ter or it has a three-dimensional one with the d z 2 charac- 



ter. These differences can be systematically understood 
in terms of the Fe-Fe and Fe-pnictogen distances. 35 ' On 
the other hand, in our calculations, the interaction pa- 
rameters do not exhibit noticeable differences between 
LaFeAsO and LaFePO. The mechanism of the supercon- 
ductivity has to clarify how these subtle differences in the 
one-body and many-body terms of the effective Hamilto- 
nian lead to a large difference in the critical temperature, 
26-55 K for As compounds and 4 K for P compounds. 

The present model may be used for further stud- 
ies by using low-energy solvers such as the dynamical 
mean-field theory 36 ' and path-integral renormalization 
group. 37 ' 38 -* Our ten-band model may also offer a firm 
starting point for further downfolding to derive effective 
models with fewer number of bands near the Fermi level. 
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